function Y = G(U)
% G - y方向通量
%   
[rho,u,v,p] = U2Variables(U);
E = U(4);
Y=[rho*v;
   rho*u*v;
   rho*v^2+p;
   v*(E+p)];